Problematic solutions to the Landau-Lifshitz 



(N 
O 

o 

cL,: equation 
< 

q\ \ W. E. Baylis and J. Huschilt 

February 2, 2008 

Oh' 

I ' Abstract 

, A critical look at the Landau-Lifshitz equation, which has been re- 

■ cently advocated as an "exact" relativistic classical equation for the mo- 

j tion of a point charge with radiation reaction, demonstrates that it gen- 

O 



erally does not conserve energy-momentum. Its failure is dramatic in the 
one- dimensional case of a stepped electric field, where it predicts discon- 
tinuous velocity and thus infinite radiation. The Lorentz-Dirac equation, 
on the other hand, in spite of its preacceleration over distances comparable 
^T*' to the classical electron radius, does not display such problems. 

1 Introduction 

(N : 

> ■ Well-known problems of the Lorentz-Dirac (LD) equation, @, |[ |, f§ |, g, g, g 

such as runaway solutions, preacceleration, nonuniqueness, have led to proposals 
of modified classical equations for the relativistic motion of a point charge in 
electromagnetic fields, including radiation reaction. Recently Rohrlichj9j has 
f^i ■ asserted that forms of the Landau-Lifshitz (LL) equation^, represent the 

Cs) \ exact equation and should replace the flawed one (LD) derived by Dirac. The 

purpose of this letter is to point out that solutions to the LL equation are not 
c/2 \ generally consistent with Maxwell's equations for the radiation. The result is a 

fairly simple to show but does not seem to have attracted much attention. 
The LD equation for the motion of a point charge e of proper velocity u in an 
J^~* external electromagnetic field F can be written in the index- free formulation Q 

based on Clifford's geometric algebra of physical space as 

mw = e((F + F se | f )u) g j , (1) 

' where the first term on the right / = e (Fit)^ = e Fu + (Fit)^ /2 is the co- 

variant Lorentz force, and F se if = mT e d{uu) /dr is identified as the effective 
field of self interaction. We use SI units with c = 1, dots indicate derivatives 
with respect to proper time, and r e is two-thirds the time for light to cross the 
classical electron radius: 

2 Kp 2 

Te = ±±^~ 6 .266 xl(r 24 s (2) 
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with K = (4tt£ ) • 

Equation (Q) can be expanded in components to give the standard tensor- 
component form m?i M = e (F^ u + F^ f ) u v , where Fj^ f = mr e d (u^u"^ / dr, 
the summation convention is adopted, the brackets [■ • ■ ] indicate the antisym- 
metric part, and the metric tensor is (i7 M „) = diag(l, —1 — 1 — 1) . However, 
the component-free algebraic formulation is cleaner and offers additional com- 
putational tools. In it, paravectors (scalars plus vectors) represent spacetime 
vectors. For example u = 7 + u = u M e M is the proper velocity with time com- 
ponent 7 = 7eo and spatial part u. An overbar indicates the Clifford conjugate 
u = 7 — u and the Lorentz-invariant square norm is uu = 7 2 — u 2 = u^u v e^e Pl 
which gives the Minkowski spacetime metric rj^ v = (e^e^g as the metric of 
paravector space, where (x) s = ^ (x + x) is the scalar part of any element x . 
Since u is a unit paravector, uu = I, and the proper acceleration u is orthogonal 
to u : (uu) s — 0. As a consequence, uu is a biparavector (a spacetime plane, 
represented by a complex vector). 

The expansion of F se if gives 

d(uu) .._ . . 

F se |f = mr e — - — = mT e uu — F (3) 
dr 

where the Lorentz-invariant Larmor power P = —mT e uu is seen from Maxwell's 
equations dF = /loj to be the power radiated by the accelerating point charge. 
It is easily seenQ that the LD equation (|l|) conserves energy and momentum 
with the radiation field between any two points on the world line of the charge 
where the acceleration u is the same. 

The LL equation^) is obtained by replacing mil in the radiation term in ([!]) 
by the Lorentz force / 



HI I'LL = I + T, ( — (fu) U 



= f + r e (f + (fu) s uj, (4) 

where the second line follows from the reality of / and its orthogonality with 
u.This is the equation given by Ford and Q'Connell[|ri), by Spohn|lC|, and by 
Rohrlich|)] in his (4a). It is easily verified that u here remains orthogonal to u. 
As Rohrlich noted, the last term in (Q) dictates that the Larmor radiation term 
P is replaced in energy-momentum conservation for the LL equation by 

ftL = -r e (/S) s . (5) 

However, this can differ from the Larmor power P, and it is P that is given by 
Maxwell's equations. 

The LL equation is the first term in an iterative expansion of the LD equation 
in powers of r e : 

mu {n+1) = 



/ + (-(*<»>«)«) (6) 
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where the nth-order approximation of u and 

^ (0) =/. (7) 

The lowest-order difference between the proper accelerations of the LD and LL 
equations is the second-order term ii — vP^ w r| (tt + uuw) . 

Rohrlich also gives an alternative form [his (4b)] , obtained by replacing / = 
e (Fu^ + e (Fu)gj with 

/ R = e( U a) s (F U ) ;R + ^(F(Fu) 3i ) SJ (8) 
The difference / — /r is first order in r e . 

2 A simple example 

As a simple example by which to compare solutions of the equations, consider 
one-dimensional motion in a pure electric field F = E = Ee, where e is a fixed 
unit vector. Express the proper velocity in terms of the rapidity w as u — 
exp (we) = cosh it; + esinhw. Then u = wue and the Lorentz force is / = 
eEu.The LD equation ([[]) becomes 

eE 

W = h T e W , (10) 

to 

whereas the LL equation (Q) has the form 

(e + t b E) . (11) 
The LL equation is the first-order iteration of the equation 



e 

TO 



to 

w& = — (12) 

TO 



(E + Ta £\ 



TO TO 
TO *■ — ' CLT K TO V CLT I 

In the limit n — > oo, the iterative solution is seen to approach the LD solu- 
tion. To lowest order, the difference between the LD and LL equations is the 
second-order term w — w^- 1 ' ~ (e/m) t^E , the power difference is Pr — P ~ 
- (e 2 /m) t^EE, and / - / R ~ (e 2 /m) T e uEE, none of which generally van- 
ishes. 
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Ford and Q'Connell]lT[ (1993) derive analytical solutions of the LL equation 
(O) for motion of a charge through an electric field in the shape of a step: 



0, x < 
E(x)={ E , < x < L 
0, L < x 



(13) 



They also show that these are the smooth limit of numerical solutions for a 
smooth rise and fall of the field. Let r = be the proper time that the charge 
enters the field from the left ( x = ) and r = T\ the proper time that it exits 
at x = L. Integration of (p"l|) gives 



W , 

W + £ + £t/t s , 
w 2 = wo + £Ti/T e 



T < 
< T < Tl 

n < t, 



where 



eE T e T e 
m L 



(14) 



(15) 



Note that wll is discontinuous: it jumps by e as the charge enters the field 
at r = and then by — e as the charge leaves at r = t\. Consequently, the 
acceleration has infinite spikes as the charge enters and leaves the field. In 
terms of Dirac delta functions S (r) and Heaviside step functions 6 (r) , 



= e [6 (t) - 6 (r - n)} + —0 (r) 9 (n - r) 

m 



(16) 



Although w and consequently m are infinite at r = 0,ti, they are integrable. 
However, the Larmor radiation, proportional to — iiit = w 2 , is not. Thus, accord- 
ing to Maxwell's equations for the field of a point charge, infinite energy is radi- 
ated from the discontinuities. The distance x traveled in the region < x < L 
is related to r by integration 



x (r) = / sinhw(r') dr' 
Jo 



L 



[cosh (wi +er/T e ) — coshwi 



(17) 



with wi — wo + £ • In particular, x (t±) = L. The energy gain of the charge is 
m (72 — 70) , where jj = cosh Wj and wi = cosh -1 (71 + a) — e. To second order 
in e, 



72 — 7o = (71 + ") cosh e 
~ a + £ 



\J (71 + a) 2 — 1 sinh e — 70 



(o) 
u — u 2 



with vffl — J (70 + a)' 



1. 



70+2 



(70 + a) up 
„.(°) 



(18) 
(19) 
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Let's compare this to solutions of the LD equation ([To|). In numerical so- 
lutions, runaways are avoided by integrating backward in time]l2"|. In analyt- 
ical solutions, they are avoided by assuming linv^oo tE (t) = and putting 
w (oo) = 0. This gives the usual integral form that includes brief periods of 
preacceleration: 

e f°° 

w (t) = — / dsE(r + T e s)e s 

( (e/rrc) E Q e T / T - (l - e^ 7 " 1 ^) , 
= I (e/m)E (1 _ e (T-n)/r^ ; 

A further integration gives 

( w + ee T ^ T - (1 -e- T ^ T ') , r<0 
io(r) = ^ w + £[l + T/r e ~e( T - Tl )/ r =] , 0<r<T! . (22) 



(20) 

r < 

< r < ri • (21) 

Ti < T 




Figure 1: Rapidity w in the case of a stepped held. Solid line: the LL solution; 
dotted line: the LD solution. Note that w is continuous in the LD solution 
by virtue of the preacceleration over times of about r e , but in the LL case, 
where there is no preacceleration, w is discontinuous. (The size of r e is greatly 
magnified for clarity.) 



3 Discussion 

Derivations of the LD equation (|l|) generally assume an expansion u in powers of 
the time difference corresponding to the effective size of the charge. The limit of 
vanishing size is then taken, traditionally with mass renormalization, although 
such renormalization can be avoided by taking specified combinations of the self 
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field [jl], ^, [lj| [l4| fl5| . One cannot expect the proper velocity to be an analytic 
function of position in regions where the field itself is discontinuous. However, 
discontinuous fields are simply idealizations convenient for finding analytic so- 
lutions. Solutions of both the LD and LL equations can be found numerically 
for more realistic field configurations, and they approach the analytic solutions 
in the appropriate limit. 

Rohrlich|| claims to have derived the LL equation as an exact classical 
equation for the point charge. However, his derivation, like most othcrsQ, 
makes the substitution of an approximate expression from the Lorentz-force 
equation with the justification that the radiation term is small. He then claims 
that because higher-order derivates of the velocity than second disappear in 
Dirac's derivation when the limit of vanishing charge radius is taken, one should 
also be able to ignore corresponding derivatives in the field. This approach 
appears to argue more forcefully for the correctness of the LD equation, which 
as seen above conflicts with the LL equation. RohrlichQ also claims that the LL 
equation has been obtained in a rigorous mathematical argument by Spohnjiof, 
but Spohn obtains his critical surface perturbatively and does not claim it to 
be exact to all orders of r e . 

The LL equation differs from the LD equation only in second order in r e and 
its solutions to realistic problems are practically indistinguishable from those 
of the LD equation since r e is orders of magnitude smaller than the smallest 
measurable time interval. Nevertheless, as seen above, it is inconsistent with 
Maxwell's equations for the radiation of a point charge and this inconsistency 
is dramatic in the case of rectiliner motion through a stepped field. This is in 
contrast to the LD equation, which is consistent. 

YaghjianQ has proposed a different "correction" to the LD equation. (Most 
of his book discusses a model of the electron as a spherical insulator of finite 
radius with a fixed surface charge, but the last section discusses the limit of 
vanishing radius to find the motion of a point charge.) He argues that the 
radiation terms do not act until the field is turned on and consequently should 
be multiplied by a scalar function that approaches a step function in the limit 
of a point charge. This eliminates preacceleration. Although his formulation 
does not explicitly treat other abrupt changes in the field, for consistency we 
assume that the sudden drop in the stepped field has no effect on the radiation 
terms until t\ when the charge leaves the field. However, this prescription when 
applied to the stepped field gives precisely the motion of the charge without 
any radiation reaction. It is therefore also inconsistent with energy-momentum 
conservation and Maxwell's equation. 

As frequently pointed out|l(|, the problems of the LD equation occur at dis- 
tance scales well below the Compton wavelength, where quantum effects become 
important. Its breakdown in the description of real particles at such scales is 
therefore not surprising. Attempts appear so far unsuccessful to find an alterna- 
tive classical equation of motion for the point charge that is free from problems 
and consistent with energy-momentum conservation and Maxwell's equations. 
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